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Abstract. We investigate the model of the disc/corona 
accretion flow around the black hole. Hot accreting advec- 
tive corona is described by the two-temperature plasma in 
pressure equilibrium with the cold disk. Corona is powered 
by accretion but it also exchanges energy with the disk 
through the radiative interaction and conduction. The 
model, parameterized by the total (i.e. disk plus corona) 
accretion rate, m and the viscosity parameter, a, uniquely 
determines the fraction of energy released in the corona 
as a function of radius and, in particular, the transition 
radius to the single-phase flow. 

Self-consistent solutions with the mass exchange be- 
tween phases display radial dependence of the parameters 
qualitatively different from the 'static' case, without the 
mass exchange. Corona covers the entire disk. The char- 
acter of the radial dependence of the fraction of energy 
dissipated in the corona is qualitatively different for low 
and high total accretion rate. 

If the total accretion rate is low, the corona becomes 
stronger towards the central object, and finally the disc 
completely evaporates, changing the accretion pattern 
into the single hot advection-dominated accretion flow 
(ADAF). For intermediate accretion rates the reverse pro- 
cess - condensation - becomes important, allowing possi- 
bly for a secondary disc rebuilding in the innermost part of 
the system. High accretion rates always prevent the tran- 
sition into ADAF, and the cold disk extends down to the 
marginally stable orbit. 

The transition radius, rtr, between the outer, two- 
phase flow and the inner, single-phase, optically thin flow, 
is equal to 4:.51m~'^^^al iRschw for to < 6.9 x 10~^Q!o;f 
and then contracts to the marginally stable orbit in a dis- 
continuous way above this critical value of to.. 

This model reproduces all characteristic luminosity 
states of accretion black hole without any additional ad 
hoc assumptions. In particular, the mechanism of the disk 
evaporation leads to a new, almost horizontal branch on 
the accretion flow's stability curve (i.e. the dependence 
of accretion rate on surface density) at the critical ac- 
cretion rate. This branch, together with the upper, ad- 
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Fig. 1. The crossection through the central region of ac- 
cretion flow into a black hole. 



vection dominated branch for optically thick disks, form 
boundaries for the time evolution of unstable, radiation 
pressure dominated disk. Therefore the disk at high ac- 
cretion rates, corresponding to Very High State in GBH 
and perhaps to Narrow Line Seyfert 1, and quasar stage 
may oscillate between the disk dominated state and the 
evaporation branch state, with only a weak contribution 
from the cold disk emission. The position of this branch 
for a = 0.1 with respect to the gas pressure dominated 
branch is consistent with the presence of only weakly vari- 
able High State in GBH and the absence of a similar state 
in AGN: all the quasars vary considerably if monitored in 
timescales of years. We also suggest a new interpretation 
of the Intermediate State, consistent with the presence of 
the strong reflected component. 

Key words: Radiative transfer, Accretion disks, Galax- 
ies:active, Galaxies:Seyfert, X-ray:stars, X-rays:galaxies 
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1. Introduction 

There are strong observational evidence for the existence 
of the hot plasma ~ lO^K in the neighborhood of the ac- 
cretion disc in active galactic nuclei (AGN) and galactic 
black holes (GBH). Soft disc photons are Compton upscat- 
tered by this optically thin plasma and form an approx- 
imately power law hard X-ray spectral component. The 
presence of the Compton reflection hump and the cold 
iron Ka line in the hard spectrum is a signature of strong 
illumination of cold disc matter by X-ray radiation (for a 
review see Mushotzky, Done & Pounds 1993, Tanaka & 
Lewin 1995, Madejski 1999). 

The location of the hot plasma is still under debate, 
but one of the attractive models postulated by Liang & 
Price (1977) and Paczyhski (1978) introduces hot dissipa- 
tive corona above accretion disc following the example of 
the solar corona, via magnetohydrodynamic waves, jets or 
other nonthermal energy transport, and then thermalized 
in corona. But such ways of energy transport occurred to 
be inefficient, only 10% of energy generated in the disc 
can be carried out towards corona. Within the frame of 
such geometry, observations seem to require that in some 
sources and/or some spectral states the innermost part 
of the disk is either disrupted and replaced by an opti- 
cally thin flow or there is a strong coronal outflow (e.g. 
Esin, McClintock & Narayan 1997, Zdziarski, Lubihski & 
Smith 1999, Beloborodov 1999). Therefore, it is important 
to determine whether any of such effects are expected on 
the basis of the physical description of the accretion flow 
(Narayan 1997). 

The disc/corona models were recently investigated fol- 
lowing two basic lines. 

The first group of models does not explain the mech- 
anism which heats up the plasma in the disk/corona part 
of the flow. The fraction of energy generated in the corona 
/cor is treated as a free parameter of the model (Haardt 
& Maraschi 1991, Svensson & Zdziarski 1994) or changes 
according to the artificial formula being inversely propor- 
tional to the distance form a central object (Esin et al. 
1997). By changing /cor and the geometrical shape of the 
hot region (continuous, patchy or clumpy corona), such 
models can be fitted to the observed X-ray spectra (for 
review see Pountanen 1999). The main conclusion from 
those models is that static corona agrees with observa- 
tions only when it does not cover the whole disc. However, 
more accurate considerations of the reflection, taking into 
account the optical depth of the Compton heated zone 
(Rozahska & Czerny 1996, Nayakshin, Kazanas & Kali- 
man 1999) or an outflow (Beloborodov 1999) may influ- 
ence this result to some extent. 

The second group of models follow Shapiro, Lightman 
& Eardley (1976) solution when hot optically thin flow ex- 
ists in the inner part of the cold disc. They assume both 
cold and hot flow to be driven via accretion and depend- 
ing on the model they take into account mechanical energy 



transport in the radial direction: like advection (Ichimaru 
1977; Narayan & Yi 1994, 1995a, 1995b; Abramowicz et 
al. 1995; Chen et al. 1995; for review see Narayan, Ma- 
hadevan & Quataert 1998 and Kato, Fukue & Mineshige 
1998), turbulent diffusion (Honma 1996), hydrodynamical 
inflow/outflow (Narayan & Yi 1994, 1995a; Blandford & 
Begelman 1999) and/or in the vertical direction: thermal 
conduction, evaporation or coronal wind (Meyer & Meyer- 
Hofmeister 1994 hereafter MMH94, Meyer-Hofmeister & 
Meyer 1999 Liu et. al 1999, TuroUa & DuUemond 2000). 

In both groups of models the crucial thing is the tran- 
sition between the hot and the cold phase. In general pic- 
ture such a transition occurs in two directions: vertical 
transition when corona accretes above a cold disc, and ra- 
dial transition when disk/corona system is replaced with 
a single optically thin flow close to the black hole (Fig. |l]) . 
Therefore, as was pointed by MMII94, in general picture 
both radial and vertical structure of the hot and cold phase 
should be considered. 

In this paper we present the simplest disc/corona 
model which is the combination of two groups of models 
described above. The similar model was presented by Esin 
et al. (1997) with better description of ADAF region, and 
with non-local radiative interaction between the hot flow 
and cold disc. Nevertheless, they were not able to deter- 
mine self consistently the transition between disc/corona 
flow, since they did not calculate its vertical structure. 
Also MMH94 stressed that thermal conduction and mass 
exchange between disc and corona in vertical direction are 
important do determine the transition. However, in their 
model (Meyer-Hofmeister & Meyer 1999 Liu et. al 1999) 
the radiative interaction between phases was neglected, an 
assumption of a single-temperature plasma was used, and 
the conductive flux on the bottom of the hot phase was 
artificially taken from the solar corona. 

Our model takes into account radiative coupling be- 
tween phases, thermal conduction and evaporation in ver- 
tical direction, and radial energy transport via advec- 
tion. The mass exchange rate is calculated self-consistently 
from the physical conditions in the transition zone. Hot 
corona is treated as two temperature plasma and, simi- 
larly to the disc, it is powered by accretion (Witt, Czerny 
& Zycki 1997; hereafter WCZ97, Esin et al. 1997). 

On each radius the vertical structure is calculated de- 
pending on comparison of accretion and vertical evapora- 
tion timescales. When accretion time scale is long enough, 
the mass exchange can be neglected and the 'static' equi- 
librium is achieved. In such a case, for vertical transi- 
tion we use the condition based on considerations of Kro- 
lik (1998), Dullemond (1999), and Rozahska & Czerny 
(2000) which is relevant to the dissipative corona. When 
accretion time scale is short mass exchange cannot be ne- 
glected, and we determine the evaporation rate through 
the disc/corona transition. 

Vertical properties are used to compute the radial 
trends solving continuity equation. The aim of this pa- 



A. Rozaiiska & B. Czerny: Disk/corona flow and transition to ADAF 



3 



per is to determine the radial dependence of the fraction 
of energy liberated in the corona for a given accretion rate, 
viscosity and a mass of the black hole and to check when 
the disk evaporation leads to a transition from disc/corona 
flow to the single optically thin flow. The radial transition 
occurs at the radius when all the mass starts to be carried 
by the corona. 

Observational consequences of our model occur in the 
dependence of the transition radius on accretion rate, 
which can be seen in different shapes of the reflected com- 
ponent (Done & Zycki 1999). Also the evaporation leads 
to the new stable branch on the accretion flow's stabil- 
ity curve ( i.e. m(S)), allowing for new evolution cycle to 
exist. If accretion rate is too low optically thick disk can 
evaporate completely in the innermost parts and the posi- 
tion of transition radius moves outward. This can explain 
the different luminosity states in AGNs and GBHs 

In Sec H we present the model. Calculation are shown 
in Sec. and we discuss our results in Sec. ^. Evaporation 
rate is determined in Appendix A. 

2. Model 

The basic problem in modeling the two-phase accretion 
flow onto a black hole is the determination of the relative 
importance of the hot optically thin flow and cold optically 
thick flow. It was shown by MMH94 in the case of dwarf 
novae that such estimates have to be based on the study 
of the vertical structure of such a flow since only in that 
case we can determine the mass exchange between the 
phases and the conditions of their coexistence. Therefore 
we follow this basic approach but we adjust the physical 
input to the conditions appropriate for the accretion close 
to the black hole. 

The main difference from a dwarf novae is that in the 
case of inner regions of the accretion disc in AGNs and 
GBH, two-temperature plasma should be taken into ac- 
count (Narayan & Yi 1995b). This is easy to see from the 
following simple argument. With the assumption of local 
viscous heating of ions, the total flux generated locally in 
the standard disk is: 



(1) 



where fix is the Keplerian angular frequency equal to 
{GM /r^Y^"^ ^ a is the viscosity, H is the height scale of 
the flow and P is the pressure. 

On the other hand, the total flux dissipated locally in 
a stationary Keplerian disk can be described as: 



^ gen — 



'SGMM 



(l-/(r/i?schw)), 



(2) 



where /(f/i?schw) represent the inner boundary condi- 
tion, in our case taken from the Newtonian approximation, 

/(r-/i?Schw) = (3i?Schw/0'/'- 



The optical depth of the flow can be described by 
— Hi^sP^} where is the opacity coefficient for elec- 
tron scattering and p is the density. We use also following 
equation of state: P = -^—pT, with the value of molecu- 
lar weight p = 0.5 for cosmic chemical composition. Boltz- 
mann constant is denoted as k and the mass of hydrogen 
atom as toh- 

Combining these equations, we can estimate the tem- 
perature of the optically thin t^s ^ 1 flow as: 

T = 2.31 X 10"!!^(i?schw/r)3/2[l - (3i?schw/r-)i/2] (3) 
a 

We use accretion rate m = M/M Edd in units of Eddington 
accretion rate described as: 



MEdd 



AnGM 



(4) 



with accretion efficiency ?/ — 1/12. Thus the temperature 
does not depend on the mass of the black hole M, and can 
be computed for assumed value of Eddington accretion 
rate, radius and viscosity. Taking a realistic value m/a ~ 
1, and the distance appropriate for dwarf nowae which 
is 10''i?schw we obtain the temperature equal to 10^ K, 
reasonable for the single temperature plasma. But, for r = 
10i?schw, the temperature attains the value of 4 x 10^^ K, 
far too high for electrons, according to the observational 
data, and two-temperature plasma should be considered. 

2.1. Basic assumptions and model parameters 

In the present paper we consider the case of stationary 
accretion onto a non-rotating black hole. We assume that 
accreting material forms an optically thick disk underlying 
an optically thin hotter skin - a corona. The flow is de- 
scribed as a continuous medium so a possible dumpiness 
of either the disk or the corona are neglected. The transi- 
tion of the flow in the vertical direction from a corona to 
a disk is gradual and corresponds to a change of the cool- 
ing mechanism with an increasing density and pressure 
towards the equatorial plane. Since the computations of 
the optically thick part of the flow required another ap- 
proach than the optically thin parts we still introduce a 
formal division of the flow into a disk and the corona but 
the physical parameters at the bottom of the corona (pres- 
sure, density and temperature) are equal to the disk values 
so there is no discontinuity between the two media. 

We assume that both the disk and the corona are pow- 
ered directly by the release of the gravitational energy 
due to accretion. Viscous torque is parameterized as in 
the standard disk of Shakura & Sunyaev (1973). This ap- 
proach is supported to some extent by the magnetohy- 
drodynamical simulations of Hawley, Gammi & Balbus 
(1996). Such assumption was used e.g. by MMII94 and 
WCZ97. Since the corona is accreting we assume the di- 
rect heating of ions and the subsequent energy transfer 
to electrons through the Coulomb coupling. Therefore, 
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the corona is basically a two-temperature medium un- 
less the Coulomb coupling is efficient enough to equalize 
the temperatures. This assumption is a simple, well de- 
termined starting point but in the future the ion-electron 
energy exchange should be considered more closely. Pro- 
ton heated plasma, adopted in our paper, is hardly ther- 
malised (Gruzinov & Quataert 1999). Spruit & Haardt 
(2000), motivated by large mean free path of ions consid- 
ered the disk heating not as conduction due to electrons, 
but conduction due to the ions. Another point is that pro- 
tons are heated efficiently only in the case of a weakly 
magnetized plasma; in a plasma with field in equipartition 
heat goes to electrons (Quataert & Gruzinov 1999). Also 
Ohmic heating stressed by Bisnovatyi-Kogan & Lovelace 
(1997) can efficiently heat up the electrons and advection 
may be not important. Diminishing the role of advection 
may have further consequences for the efficiency of mass 
exchange between the disk and the corona (Rozanska & 
Czerny 2000). 

We take into account the radial advection, Compton 
and bremsstrahlung cooling and the conduction due to 
electrons and the vertical mass exchange, which are im- 
portant at the basis of the corona. Disk and corona are 
radiatively coupled, as in the classical paper of Haart 
& Maraschi (1991) although the cooling process is now 
treated locally. We assume that, the cold disc is the origin 
of the seed photons at a given radius. We neglect the con- 
tribution of the synchrotron radiation to the seed photons 
for Comptonization (Narayan & Yi 1994), as well as non- 
local photons which were included in the study of Esin et 
al. (1997). We also neglect the formation of wind at the 
top of corona discussed in some papers (e.g. Begelman, 
McKee & Shields 1983, Begelman & Blandford 1999). 

In flows cooled by advection, the motion may be not 
Keplerian and proper calculation of equation of motion 
should be done to describe this effect. Nevertheless, the 
aim of this paper is to calculate properly the vertical struc- 
ture while the radial one is computed less carefully. There- 
fore, together with hydrostatic equilibrium in vertical di- 
rection, we assume the Keplerian motion of the material. 

The model parameters are: mass of the black hole, M, 
accretion rate, to, and the viscosity parameter a, assumed 
here for simplicity to be the same both in the disk and 
in the corona. All other quantities, including their radial 
dependences, are uniquely determined by those global pa- 
rameters for a stationary flow. Corona model does not 
actually depend on M but there is a weak dependence of 
the cold disk properties on the mass of the black hole itself 
so this parameter cannot be entirely neglected. 

2. 2. Vertical structure of the corona at a given radius 

In the present study we simplify the problem of hydro- 
static equilibrium. Instead of solving complex differential 
equations as in MMH94, we simply study algebraic de- 
pendences of various quantities on the pressure which is a 



monotonic function of the distance from equatorial plane. 
With such assumption the local viscous heating of ions 
can be described following the Eq. |l|: 

Q+ = ^aPQK, (5) 

A fraction of this energy is carried by advection. This 
fraction is proportional to the ratio of the ion tempera- 
ture Tj to the virial temperature with the coefficient S 
determined by the radial derivatives of thermodynamical 
quantities (e.g. Ichimaru 1977, Muchotrzeb & Paczyhski 
1982, Abramowicz et al. 1995, Janiuk, Zycki & Czerny 
2000). At present we simply assume 6=1 but in principle 
this value should be determined iteratively. 

Remaining fraction of energy is transfered to electrons 
via Coulomb coupling 

Q+ f 1 - 5^) = Qt = ^^—Dp^m T^)T-'/' (6) 
V Tv y 2 TOH 

where is the electron temperature, and D is the 
Coulomb coupling constant equal 4.88 x 10^^ (in s^^). At 
this exploratory stage, we neglect relativistic corrections 
to this equation given by Stepney & Guilbert (1983). 

In the upper layers of the corona the electrons cool 
down radiatively, i.e. 

Q+=A(p,Te). (7) 

The processes included are the Compton cooling due to 
the presence of soft photons -Fgoft from the disk and 
through the bremsstrahlung cooling 

A(p, Te) = FsoftKcs^pTe + Bp^Tl'^. (8) 

Here -Fgoft is the soft radiation flux coming from the disk 
due to the energy generation in the optically thick part of 
the flow as well as due to interception and reemission of 
the coronal emission Fhaid (Haardt & Maraschi 1991) 

^^soft = Fdisk + 0.5Fhard, (9) 

TOe - mass of the electron, c - light velocity, and B 
- bremsstrahlung cooling constant equal 6.6 x 10^° erg 
s^^cm^^ g~^. We assume that half of the coronal flux 
Fhard, illuminates the disc and all of this is absorbed and 
subsequently reemited, so albedo is equal to zero. 

In Eq. 1^ we include only the Compton cooling and we 
neglect, for simplicity, the Compton heating term which 
may be important within the transition zone where the 
electron temperature drops below 10* K. 

Equations |^ and ^ determine the temperature profile 
in the upper part of the corona, or more precisely, the 
dependence of the ion and electron temperature on the 
gas pressure which we identify with the ion pressure for 
simplicity, P = k/mnpTi. 

At some depth, however, a transition to a cool disk 
is expected. Within such a transition zone the essential 
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role is played by conduction, neglected so far. Here we as- 
sume that the transition from hot to cold gas happens very 
sharply, under the constant pressure. Phase, so the drop in 
the temperature is simply compensated by the increase in 
the density. Determination of this pressure is the key ele- 
ment of the model since the value of this pressure uniquely 
determines the fraction of the energy /cor liberated in the 
corona in the following way. 

Since the ion temperature is close to the virial tem- 
perature, the pressure scale high is roughly of the order 
of the radius, r, so instead of solving precisely the hydro- 
static equilibrium we can estimate the total energy flux 
generated in the corona as 



Pcor Q[ (Pbasc)^, 

and the flux emitted by corona as 



(10) 



(11) 



The total energy dissipated due to accretion both in the 
disk and the corona is given by formula: 



disk 



0.5Pha 



1 



6^ 



1 



5^ 



(12) 



which relates to the accretion rate both within the disk 
and the corona due to the assumption of stationarity and 
Keplerian motion and Ftot = Pgon given by Eq. 

In the accreting corona model, the fraction of energy 
dissipated in the corona corresponds to the same fraction 
of the mass accreting through the hot phase, i.e. 



(13) 



In this way the determination of the pressure at the 
basis of the corona allows to determine its structure at 
a given radius from the global model parameters, M, m, 
and a. Therefore the way how we determine this value is 
a crucial element of the model itself which removes arbi- 
trariness from the disk/corona accretion flow present in 
most models apart from MMH94, Zycki, CoUin-Souffrin 
& Czerny (1995), WCZ97 and the directly related papers. 

The computations of the vertical structure of the op- 
tically thick part of the flow are not necessary for the 
computation of the corona. However, when the corona 
structure is determined the disk structure is calculated, 
with the boundary conditions imposed by the presence of 
corona, as described by Rozahska et al. (1999). No itera- 
tions between the disk and the corona computations are 
needed. 

The properties of the transition zone are calculated fol- 
lowing the approach of Krolik (1998), DuUemond (1999), 
and Rozahska & Czerny (2000). We consider two cases: 
a 'static' balance solution without the mass exchange be- 
tween the disk and the corona and the self-consistent case 
based on the continuity equation of the coronal flow. 



2.2.1. Equilibrium solutions without mass exchange 

Interesting and important solutions are obtained under an 
assumption that the accretion flow is slow in comparison 
with the timescale of exchange of the mass between the 
disk and the corona. In that case the equilibrium solution 
is defined by the condition that there is no mass exchange 
between the disk and the corona (see Rozahska & Czerny 
2000). 

To describe the equilibrium solution we integrate the 
energy balance equation as shown in Appendix A, with 
velocity in vertical direction equal zero (see Eq. A. 3). 
It meant, that the heat flux which is transported via elec- 
trons from hot upper corona towards the cold matter is 
balanced by radiative cooling operating in the disc. Since 
the transition is described under constant pressure the 
appropriate criterion reduces to an integral (Rozahska & 
Czerny 2000) 



[Qt - K)KoT^'^dT, = 



(14) 



where Tg is the cold disk temperature and T'e(P) is the 
electron temperature within the upper corona part deter- 
mined at the pressure P and satisfying the heating/cooling 
balance formulated by Equation 0. Here, n — KqT^/^ is 
thermal conductivity coefficient in ergs cm~^s~^K^^ and 
conductivity constant Kq is assumed to be equal 5.6 x 10^ ''. 
The ion temperature in this equation has to be determined 
from the Coulomb coupling equation (Eq. ^) and the den- 
sity from the value of the (constant) pressure and the ion 
temperature. 

Equation |l^ determines the pressure of the transition 
zone, thus closing up the set of equations determining the 
disk/corona flow. 

2.2.2. Solutions with evaporation/condensation 

The mass exchange between the disk and the corona can- 
not be neglected if the timescale of accretion is short so 
the 'static' type of equilibrium cannot develop. It is also 
necessary if we want to describe the process of corona for- 
mation. This approach is more appropriate but more com- 
plex since it is essentially non-local and requires a global 
solution for all radii. Here we basically follow an approach 
of MMH94 by finding a family of numerical solutions at a 
single radii, determining their scaling properties with the 
radius and finally finding unique global solution at the 
basis of this scaling. 

Therefore we allow the Equation |l^ to be not satis- 
fied which means that we have to solve the vertical en- 
ergy transfer in the transition zone allowing for evapo- 
ration or condensation, depending on the choice of the 
pressure Pbaso- If radiative processes are to low to balance 
the heat flux, the disc will evaporate as shown by DuUe- 
mond (1999). However, radiative cooling may overwhelm 
the heat flux if pressure at the base is high enough (higher 
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than the pressure under which mass exchange equals to 
zero, as shown by McKee & Begelman in 1990), and in 
such situation matter will condense. 

To determine the mass exchange rate we integrate the 
energy balance equation ( A.l) in the way detaily presented 
in Appendix A. Combining equations (A. 4) and (A. 7) the 
evaporation rate rhz across the disk/corona boundary is 
given by the formula. 



TOz 



A] KoTc^^dTc 



2 rrifi 



(15) 



The sign of this quantity is either positive (disk material 
evaporates into the corona) or negative (coronal material 
condensates into the disk), depending on the net effect of 
the electron heating and cooling balance. 

The radial dependence of Phase is subsequently deter- 



mined from the continuity equation (see Section 2.3.2) 



2.3. Radial component of the continuity equation and 
global solutions 

2.3.1. Solutions without mass exchange 

In this case the equilibrium is built locally so the deter- 
mination of the disk/corona structure at a given radius is 
completed if the local solution at this radius is found. The 
accretion rate within the corona is determined at each ra- 
dius separately and the mass exchange between the disk 
and the corona account for its dependence on the disk ra- 
dius, i.e. its evolution with the flow. All parameters, like 
accretion rate in the corona, soft X-ray flux from the disk, 
the ion and electron temperature profile, the density pro- 
file and the corona optical depth as well as the pressure at 
the basis of the corona where the temperature drops down 
to match the disk values are determined by the choice of 
the radius r and the three global model parameters, M, 
rh and a. 

2.3.2. Solutions with evaporation/condensation 

If the timescale of building the evaporation/condensation 
equilibrium is long in comparison with the viscous 
timescale in the corona we have to describe the build-up 
process of the corona, i.e. the accretion rate in the corona 
at any radius r will refiect the amount of material which 
evaporated from the disk at all radii larger and equal to 
r. 

The change of the accretion rate in the corona A/cor is 
therefore due to the mass exchange with the disk 



dM„ 



dr 



= — 27rrTO, 



which in dimensionless units reduces to 
drhcox _ 2r]KosG 



dR 



-MRrh^ = -27.9MBmRm^ 



(16) 



(17) 



where R = r/i?schw and Mbhs = M/IO^Mq. 

In order to solve this equation we follow the approach 
of MMH94. We find a family of solutions at a given radii 
for a set of values of the pressure Phase- We express those 
solutions through the computed rhcoi, thus obtaining a 
numerical relation mz(mcor)- We study the scaling of this 
relation with radius, thus finding an analytical expression 
valid in the entire disk. This analytical expression is sup- 
plied to Equation |l^ and the final analytical result for the 
dependence mcor('") is found. 

2.4- Cold disk structure 

Determination of the corona properties, i.e. the fraction 
of energy dissipated in the corona and the pressure at the 
corona basis provides the boundary conditions necessary 
to calculate the vertical structure of the cold disk flow 
(e.g. Rozahska et al. 1999). Therefore, the description of 
the hot flow influences the properties of the cold part of 
the flow to their strong coupling through these boundary 
conditions. 

In the present paper we concentrate mostly on the hot 
flow properties and the cold disk properties are only re- 
quired for the qualitative discussion of the connection be- 
tween the single phase hot flow and disk / corona flow (see 
Section 3.3) so we simplify the description of the cold flow. 
We apply the vertically averaged approach but we deter- 
mine the dimensionless coefficients involved in this pro- 
cess on the basis of full computations of the disk vertical 
structure, as in Janiuk et al. (2000) (see also Muchotrzeb 
& Paczynski 1982, Abramowicz et al. 1988). 

We assume the same value of the viscosity coefficient 
in the hot and in the cold flow although they may in 
principle be considerably different due to the different 
plasma properties. We include advection in our descrip- 
tion of the flow, and the appropriate coefficients used are: 
Bi = B2 = B3 = B4 — 1, and q^dv = 2. The disk is 
assumed to be optically thick so the current approach 
does not allow yet to follow the final transition from the 
disk/corona fiow to a single phase optically thin accretion. 

3. Results 

3.1. Solutions without mass exchange 

3.1.1. Vertical structure of the flow at a given radius 

We first discuss the vertical ion and electron tempera- 
ture and density profiles at a representative radius r = 
10.8Pschw for mass of the black hole equal lO^M© and 
the viscosity parameter a = 0.1. The virial temperature 
corresponding to this radius is equal 10^^ K. 

The set of algebraic equations (Eq. || and Eqs. ^ - ^2|) 
describing the upper layers of the corona formally posses 
two solutions for a given pressure (corresponding to a cer- 
tain distance z from the equatorial plane) : one solution is 
mostly advection-dominated, low density branch and the 
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other one is high density, mostly radiatively cooled branch. 
This property, found also by DuUemond (1999), is a direct 
analogy of the vertically averaged optically thin solutions 
found by Abramowicz et al. (1995). Both solutions exists 
for low values of the pressure but they approach each other 
and finally merge at a certain value of the pressure, corre- 
sponding to an ion temperature T; = 0.5Tv. Dipper inside 
there are no solutions, as in the case of too high accretion 
rate in optically thin flows. 

However, the situation of the disk/corona accretion is 
not the same as in the case of a single phase flow. Only 
one of the above solutions is acceptable. At the advection- 
dominated branch the ion temperature decreases with the 
increase towards the midplane and the density increase 
thus opening a possibility to match this solution to the 
cool disk part of the flow. The other branch shows opposite 
trends and therefore we reject it as unphysical in the case 
of two-phase accretion. 

The fact, that the algebraic solution does not extends 
very deeply does not pose a problem since we have now to 
take into account the transition to the disk flow under no 
evaporation conditions. The value of the pressure which al- 
lows the coexistence of the disk and coronal flow under no 
evaporation/condensation condition given by Equation ^ 
is always smaller than the value of the pressure at which 
the two solutions merge. 

The example of the change of the ion and electron 
temperature with a pressure increasing towards the equa- 
torial plane is shown in Fig. ^ for the accretion rate 
m = 1.3 X 10^'^ in units of Eddington accretion rate. 
The ion temperature is very close to the virial value, most 
of the heat generated in the corona is advected towards 
smaller radii. The electron temperature is always signifi- 
cantly lower than the ion temperature since the Coulomb 
coupling is nor very efficient but it slowly increases to- 
wards the equatorial plane (i.e. increasing P). For the as- 
sumed parameters the equilibrium pressure Pbaso is equal 
7.94 X 10^ g cm-3. 

The fraction of energy generated in the corona is equal 
0.89, so the soft photon flux is predominantly due to the 
fraction of corona emission absorbed by the disk (i^soft = 
5.2 X lO^^erg s~^cm~^). The optical depth of the corona 
is very low, 6 x 10"'^ and the density at the upper edge of 
the transition zone is equal 5 x 10~^^ g s~^ cm~^. 

The transition zone itself is also complex. We assumed 
that the transition happens at a constant pressure so its 
structure cannot be shown in Fig. ^. Instead, we draw 
the relation between the ion and electron temperature 
within this zone, as predicted, neglecting thermal conduc- 
tion (only for radiative processes and advection), in Fig. ||. 

In the upper part of the transition zone itself the ion 
temperature is still very close to the virial temperature, 
with advection as a dominant cooling mechanism for ions. 
However, deeper inside the ion temperature decreases and 
the density increases. Radiative cooling becomes more ef- 
ficient. The advection-dominated flow is replaced by the 
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Fig. 2. The change of the ion and electron temperature 
with pressure increasing towards equatorial plane for rfi = 
1.3 X 10"'^. Vertical continuous line mark the Phase value 
at which the temperature drops dramatically to match the 
disk values without, however, mass exchange between the 
disk and the corona. In such a solution /cor = 0.89. If anal- 
ogous transition takes place at lower pressure it would lead 
to continuous disk evaporation due to conduction while 
the delayed transition at higher pressure would result in 
continuous corona condensation into a cold material and 
in each case /cor would take different value. 

radiatively cooled flow, basically of the type studied by 
Shapiro, Lightman & Eardley (1976). This branch was 
found to be unstable (Piran 1978, Narayan & Yi 1995b) 
and it does not provide us with a monotonic decrease 
of the electron temperature so we assume that actually 
a sharp transition takes place along a vertical line, as 
marked in Fig. |^. Our integral of the difference between 
the heating and cooling multiplied by the conduction flux 
(Equation |l^ is calculated applying such a profile. Full so- 
lution of differential equations with the conduction due to 
electrons as well as ions would produce a smooth solution 
matching the upper and lower branch without the third, 
intermediate branch. We intend to study this problem in 
the future but for the determination of the basis solution 
properties the detailed knowledge of the transition zone is 
not necessary as the optical depth of this zone is generally 
very small. 

The fraction of the energy generated in the corona for 
an accretion rate of to = 1.47 x 10"'^ is high but generally 
this is not the case. In Fig. ^ we show the dependence of 
this value on the accretion rate at the radius 10 i?schw For 
large accretion rates the corona is weak and its strength 
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Fig. 3. The relation between the ion and electron temper- 
ature within transition zone. 



increases rapidly when the accretion rates drops down to 
m = 1.53 X 10^^. Below 1.47 x 10^^ there are no solu- 
tions for disk/corona equilibrium without the mass ex- 
change. Solutions with negative value of /disk are unphys- 
ical. Those results arc qualitatively very similar to the 
results obtained by WCZ97 which were based on approx- 
imate description of the disk/corona equilibrium through 
the condition on the value of the ionization parameter S 
(shown in Fig. ^ as a dashed line) . 

Since the corona is not isothermal there is no single 
value of the ion and electron temperature for a given ra- 
dius. However, the representative values are those at the 
top of the transition zone, at Pbaso- We plot them also 
as functions of the accretion rate (Figs. |). We see that 
the trends are similar to those determined by WCZ97 but 
there are quantitative differences. The ion temperature is 
always very close to the virial temperature and, at that 
radius, it never approaches the electron temperature. Fi- 
nally, we plot the optical depth of the corona at the same 
radius (Fig. |). The corona is much thinner than predicted 
by WCZ97, its optical depth strongly depends on the ac- 
cretion rate. 



3.1.2. Radial trends 

All quantities depend not only on the accretion rate but 
also on the disk radius. Therefore we study now the radial 
dependences of the coronal parameters for a representative 
value of the accretion rate m — 0.05. 

There is an essential, qualitative change with respect 
to the results obtained by WCZ97. 



Fig. 4. The dependence of the fraction of energy gener- 
ated in the disk /disk = 1 — /cor on the total accretion 
rate at lOiischw^ equilibrium solution (continuous line with 
squares), solution with evaporation/condensation (long 
dashed line) and prediction from WCZ97 (short dashed 
line), Solutions below dotted line are unphysical, but they 
are shown to demostrate the mathematical reason for the 
disappearance of solutions for low m. 



In WCZ97 the strength of the corona increases out- 
wards: corona is weak close the black hole while all energy 
is generated within the corona at a certain critical dis- 
tance Tmax- Beyond this radius the disk/corona configura- 
tion does not exists. According to our new description of 
the disk/corona transition under no evaporation condition 
the solutions at radii larger than a certain rmax also do not 
exists but at a radius only marginally smaller than Tmax 
the fraction of energy generated in the corona is not 1, 
but about 0.5 (Fig. In a very narrow radial region be- 
low Tinax two solutions cxist: one with weaker corona and 
the other one with strong corona, which merge at Tmax- 
The strong corona solution reaches the value /cor = 1 at 



a radius not much lower than 



and the solution does 



not continue down to smaller radii (mathematically, it ex- 
ists as unphysical solution with /cor > !)■ Therefore this 
branch of solution does not offer a global solution for the 
corona. On the other hand the weak corona branch contin- 
ues down to the marginally stable orbit and this branch 
represents the solution acceptable under adopted condi- 
tions (no evaporation/condensation). The nature of this 
multiple branch solutions is connected with the presence 
of advection, as studied for the case of WCZ97 condition of 
disk/corona transition by Janiuk, Zycki & Czerny (2000). 
The advection dominated branch tends to vanish if the 
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Fig. 5. The dependence of the ion and electron tempera- 
ture at the top of the transition zone on the total accretion 
rate at 10i?schw: equilibrium solution (continuous line with 
squares), solution with evaporation/condensation (long 
dashed line) and prediction from WCZ97 (short dashed 
line). 
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Fig. 7. The dependence of the fraction of energy gener- 
ated in the corona on the disk radius for the total accretion 
rate m = 0.05: equilibrium solution (continuous line with 
squares), solution with evaporation/condensation (long 
dashed line) and prediction from WCZ97 (short dashed 
line). 
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Fig. 6. The dependence of the optical depth of the corona 
on the total accretion rate at 10i?schw: equilibrium solu- 
tion (continuous line with squares), solution with evapora- 
tion/condensation (long dashed line) and prediction from 
WCZ97 (short dashed line). 



Fig. 8. The dependence of the fraction of energy generated 
in the corona on the disk radius for the total accretion rate 
m ~ 0.05 - version of Figure ^ expanded near the outer 
edge of the coronal equilibrium solution. 
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Fig. 9. The dependence of ion and electon temperature 
at the top of the transition zone on the disk radius for 
the total accretion rate m = 0.05: equihbrium solution 
(continuous line with squares), solution with evapora- 
tion/condensation (long dashed line) and prediction from 
WCZ97 (short dashed line). 

solution is computed more accurately by determination of 
the coefficient S iteratively from the radial derivatives of 
physical quantities. 

Therefore, according to our model, the corona in equi- 
librium with the disk covers only a finite part of the 
disk, its strength increases outwards but it never dissi- 
pates more than about a half of the gravitational energy 
available at a given radius (Fig. |^). The transition from 
outer part without a coronal solution to the inner part 
with corona is sharp if studied within the frame of so- 
lutions with no disk/corona mass exchange and the true 
nature of this transition can be understood only within the 



frame of the solutions with evaporation (see Section 3.2) 



The ion temperature at the top of the transition zone 
is always quite close to the local virial temperature so 
it clearly decreases with the radius (Fig. ||). On the other 
hand the electron temperature at Pbaso increases outwards 
on the basic branch, as in WCZ97. Its value at the outer 
edge of the coronal solution is of the same order as pre- 
dicted by WCZ97 but closer to the central object is much 
cooler. 

The optical depth of the corona is not constant, as in 
WCZ97, but it strongly decreases inwards along the basic 
branch (Fig. [Io| ). Therefore, the effect of Comptonization 
apart from the very outer edge of the coronal solutions is 
completely negligible as long as the inner disk is not evap- 
orated (see Section |3^ ). In this paper we neglect Compton 



Fig. 10. The dependence of the optical depth of the 
corona on the disk radius for the total accretion rate rh = 
0.05: equilibrium solution (continuous line with squares), 
solution with evaporation/condensation (long dashed line) 
and prediction from WCZ97 (short dashed line). 



heating, since the accretion is a main source of heating. 
However, Compton heating may increase the optical depth 
of the base of the corona and we plan to study this effect 
in a future work. 

The radial trends change with the accretion rate is a 
way similar to that described by WCZ97. The extension of 
the corona, rmax, is large for large accretion rate and de- 
creases with decreasing accretion rate. For accretion rates 
below rh = 1.68 x 10~^ no equilibrium solutions are ob- 
tained. Again, the nature of this phenomenon can be un- 
derstood by analyzing the solutions which allow for mass 



exchange in the disk/corona system (see Section 3.2). 



3.2. Solutions with evaporation/condensation 

Construction of those solutions is more complex as the 
solutions have global character, i.e. the results at a given 
radius are coupled to those at other radii through the con- 
tinuity equation for the corona (Equation |l^). Since we 
follow semi-analytical approach to solution of Eq. [l^ we 
first present local numerical solutions which are parame- 
terized by the coronal accretion rate, TOcod in addition to 
global parameters M, m and a, and we give the analyti- 
cal formulae approximating the relation between rhcor and 

The vertical structure of the disk/corona system with 
evaporation is qualitatively similar to the case without 
mass exchange. However, the global trends are essentially 



A. Rozaiiska & B. Czerny: Disk/corona flow and transition to ADAF 



11 



L5xl0-f 



C\2 

a 



10-3 _ 



.g 5x10-'° - 




0.2 0.4 0.1 

™cor/ ™ 



2x10-15 



lO-i'i - 



-IxlO-is - 



-2xl0-"i 




2x10-8 4x10- 

™cor/ ™ 



6x10-0 



Fig. 11. The local relation between the accretion rate in 
the corona and the evaporation rate of the disk at 10i?schw 
for the total accretion rate m = 5 x 10~^. Dashed line 
shows the predictions from analytical formula ^ 

different so in this section we concentrate on those aspects 
of the model. 

In Fig. ^ we show the dependence of the evaporation 
rate of the cool disk on the accretion rate in the corona 
at 10i?schw for a small value of the total accretion rate 
(m = 5 X 10^^). From Fig. ^ we see that we should not 
expect an equilibrium solution without evaporation so it 
is not surprising that the evaporation rate rhz is positive 
for all values of the accretion rate within the corona from 
to the total accretion flux. 

The relation is well approximated by an analytical for- 
mula 

= 1.27 X 10"2m^/i^i?"3/2^o;y^MBj^g [gs^^cm^^j^ (18) 

Here, viscosity is described in dimensionless units: ao.i = 
a/ 0.1. This formula neglects the effect of condensation and 
it is a good approximation of the solution only for very 
low accretion rates to. Applicability of this asymptotic 
expression will be qualitatively discussed at the end of 
this Section. 

Such an analytical expression can be now conveniently 
used to determine the global solution by substituting it to 
the continuity equation 

The radial dependence of the coronal accretion rate 
and evaporation rate are therefore: 



;.23 X lO-^MBj^sa^Y^^""^* 



(19) 
(20) 



Fig. 12. The local relation between the accretion rate in 
the corona and the evaporation rate of the disk at 10i?schw 
for the total accretion rate to = 0.05. Short dashed line 



shows the predictions from analytical formula 18 and long 



dashed line from 22 



These formulae are qualitatively similar to the formu- 
lae derived by MMH94 in the case of a single-temperature 
plasma and no condensation condition at the basis of the 
corona. The fraction of mass carried by the corona in- 
creases inwards. Finally, at a certain radius i?cvap all the 
mass is in the coronal flow due to the complete disk evap- 
oration: to — TOcor (Fig- long dashed line). Its value 



4.51ao.i^TO 



-4/3 



(21) 



is surprisingly similar to the value determined by Liu et 
al. (1999) on the basis of Mayer & Mayer-Hoffmeister's 
results. The only essential difference is the dependence on 
the viscosity present in our approach but absent in theirs. 
This is basically the result of our different description of 
the conduction flux close to the disk surface and the pres- 
ence of the additional step in energy transfer in the form of 
ion-electron coupling which results in the electron temper- 
ature much lower than the virial temperature, depending 
also on the viscosity. 

Assuming the viscosity parameter a = 0.1 we obtain 
practically the same results as Liu et al. (1999) but larger 
(smaller) viscosity in our case result in larger (smaller) 
purely 'coronal' region. The dependence on the accre- 
tion rate is almost identical in both cases (see also Sec- 
^). 



tion 



We now proceed to the case of large accretion rate. 
From Fig. ^ we see that for large accretion rate to = 0.05 
at 10 i?schw there is an 'static' equilibrium solution so 
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Fig. 13. The global solution for the radial dependence of 
the fraction of energy generated in the corona for the total 
accretion rate m = 0.05 and a = 0.08 (thick continuous 
line). Short dashed line shows the predictions from ana- 
lytical formula ([T9|), long dashed line from an analytical 
approximation to numerical results shown in Fig. ^ as a 
continuous line with squares, and dotted line shows the 
solution of WCZ97. 



generally we can expect either evaporation, or condensa- 
tion, depending on the value of the accretion rate within 
the corona itself. Fig. |l^ confirms those expectations. For 
extremely low values of coronal accretion rate rhz is still 
well approximated by the analytical formula (Eq. ^8|) but 
a turnover appears at larger values and finally evaporation 
is replaced with condensation, as becomes negative. 

Therefore we need to replace Equation |l^ with more 
general formula able to represent both evaporation and 
condensation phenomenon. 



1.27 X 10-2m5/3i?-3/2ao.i-7/2M^^^ X [1 



-1 X 10^^m^°Th''^^^R 



cor 
5/2 n-12.5 



ao.i 



-8.5 



BH8 

[gs^-^cm^^] 



(22) 



This formula is oversimplified as it does not take into 
account the non-monotonic two-value dependence of /cor 
on the radius for models in equilibrium, for which rhz = 
(see Fig. ^). However, it reproduces accurately enough the 
extension of the strong condensation zone. More accurate 
description would require numerical computations of the 
vertical structure of the corona together with integration 
of Eq. 

Previously found global solution based solely on evap- 
oration applies everywhere in the disk down to the 
marginally stable orbit if the second term in Eq. B3 is 



smaller than 1, i.e. when the total accretion rate is below 
the limiting value rh = 1.31 x 10~'^Q;gY^^ for i? = 3. In 
objects with such a low accretion rate the corona forma- 
tion proceeds monotonically, corona carry an increasing 
fraction of mass as the flow approaches a black hole and 
finally, at radius given by Eq. |2l| the disk disappears, as 
it is replaced by an optically thin flow of ADAF type. 

In objects with significantly higher accretion rate the 
situation is more complex than in low accretion rate case. 
The equations |2| and |l^ have to be now solved numer- 
ically. An example is shown in Fig. ^ (thick continuous 
line). At the outermost part of the disk the corona forma- 
tion proceeds according to Equation |l^, as before. Closer 
in the cooling becomes more efficient and the corona for- 
mation proceeds more slowly than predicted by Equa- 
tion |l^. When we approach the radius for which an equi- 
librium solution exists (i.e. possible balance between the 
disk and the corona without evaporation, see Sect. 3.1) 



the corona reaches its maximum strength. In the inner- 
most part of the disk the corona strength decreases as the 
evaporation is replaced by condensation. The radial de- 
pendence of the fraction of energy generated in the corona 
in this region is qualitatively similar to the result given un- 
der assumption of no mass exchange between the disk and 
the corona. Trends are the same as those given by WCZ97 
although the actual radial dependence is much steeper. 

A sequence of solutions showing the radial dependence 
of the fraction of the energy generated in the corona for 
several values of the accretion rate is shown in Fig. |lj. For 
significantly higher accretion rates the maximum strength 
of corona never approaches the value of 1, so the two- 
phase flow extends to the marginally stable orbit. When 
the total accretion rate is intermediate the evaporation is 
efhcient enough to cause a transition to an optically thin 
flow (/cor = 1) at a radius approximately given by Eq. pl| . 
However, in the innermost region the condensation is effi- 
cient if there is large soft photon flux available to cool elec- 
trons through Comptonization. We expect that it might 
lead to cool clump formation or even a reconstruction of 
the disk in this region. Such a secondary disk rebuilding is 
now present in our model for accretion rates higher than 
m""". More accurate numerical results for an accretion 
rate which leads to an ADAF type flow in the innermost 
part give the value 



1.9 X lO^-'a, 



3^17/25 



0.1 



(23) 



However, the current version of our model does not 
allow to describe yet the details of the transition from 
optically thick disk through optically thin disk to complete 
disappearance of the disk flow and eventual reverse of this 
process in the innermost part of the flow. Therefore, at 
present, we usually assume that once the transition occurs 
from the flow from disk/corona geometry to an optically 
thin flow at a certain radius it means that there is no disk 
flow closer to a black h ole. However, we also discuss an 
opposite case in Section 
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Fig. 14. The global solution for the radial dependence of 
the fraction of energy generated in the corona for a = 0.05 
and the accretion rates equal 0.5, 0.05, 6.3 x 10^"^, 3 x 10^'^, 
1.13 X 10-^, 3 X lO"'', and 5 x 10"^ Dashed line shows 
the formal continuation of the solution for accretion rate 
3 X 10"'^ - reconstruction of the disk due to condensation. 



Fig. 15. The dependence of the transition radius from 
disk/corona flow to an ADAF flow on the accretion rate 
for the viscosity parameter a = 0.1 (continuous line), 
a = 0.08 (dashed line) and a — 0.05 (long, dashed line). 
Dotted lines show the corresponding predictions based on 
the pure evaporation formula (Eq. ^l|). 



Further decrease of the accretion rate results in a 
smooth expansion of the ADAF zone. In all other mod- 
els the first appearance of the single phase hot flow is ex- 
pected at the marginally stable orbit (e.g. Esin et al. 1997). 
In our model the disk either extends to the marginally sta- 
ble orbit or down to a radius larger than the minimum ra- 
dius dependent on the viscosity, with no solutions for the 
transition radius ever expected between these two radii. 
This property of our solutions is caused by the domina- 
tion of the condensation over evaporation in the innermost 
part of the disk where the pressure is large. 

The dependence of the transition radius on the accre- 
tion rate is shown in Fig. |l5|. We see that for low accretion 
rates the position of the transition radius depends on the 
accretion rate and the relation is well described by Eq. 
as expected. However, at a certain value of the accretion 
rate the condensation takes over and the coronal dissipa- 
tion decreases in the innermost part of the disk enough 
to prevent the transition to an ADAF flow. This critical 
accretion rate Wcvap. branch is well described by the two an- 
alytical conditions: Eq. |2l| and rhz = (see Eq. |2^) which 
leads to 



Wcvap. branch = 6.92 X 10 ^ao^l- 



(24) 



This accretion marks the formation of the evaporation 
branch of the disk/corona solution which is essential for 
explanation of specific luminosity states of accreting black 



holes and time-dependent behaviour of X-ray sources (see 
also Sec. 



3.3) 



However, the essential point which directly results 
from radial computations is the understanding of the in- 
ner region without disk/corona solutions. What WCZ97 
interpreted as a 'bare' disk is actually a 'bare corona', or 
ADAF type of flow. 

Determination of the coronal accretion rate as a func- 
tion of the total accretion rate is a starting point of the 
discussion of the stability of solution and eventual time- 
dependent evolution. 

3.S. Expected time- dependent behaviour 

3.3.1. Surface density of the accretion flow 

The stability and time evolution of the accretion flow is 
most conveniently discussed on the plot showing the ac- 
cretion rate, m, and the surface density of the flow, E, at 
a representative radius. 

We calculate the fraction of the energy generated in the 
corona at by numerical integration of Eq. |l^, assuming the 
analytical formula for the local evaporation/condensation 
rate given by Eq. Knowing the corona properties we 
solve the cold disk structure at this radius, with the ap- 
propriate boundary conditions (see Sec. 2.4). 
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Generally, for rh > rrievap. branch the strength of corona 
is small and the solution is practically the same like in the 
thin disc without corona. However, when rh approaches 
Wevap. branch this is no lougcr true and the solution starts 
to depart horizontally from standard log m-log S curve 
(Fig. 16). The plot depends on a distance from the black 
hole, i.e. whether it is higher or lower than the critical 
radius for which the curve in Fig. |l| has maximum at 
/cor = 1- The critical radius can be determined analyti- 
cally by combining Eq. 21 and Eq. |2j which gives 



(25) 



For R > i?cvac, /cor incrcascs with decreasing m, and 
finally /cor = 1 for a certain value of accretion rate rhcvap 
smaller then mcvap. branch- As /cor ^ 1, the disc surface 
density drops with m —> movap- 

For R < Rcvac the plot depends on the assumption 
of considering or rejecting the solutions representing the 
secondar y d isc condensation (see dashed line in Fig. ^ 
and Sec. |3.2|) . In Fig. |l^ we show the relation between the 
accretion rate and the surface density for r = 10i?schw and 
a = 0.1. 

If we assume that once the disk is evaporated it never 
forms again closer to the black hole, then the transition 
radius between the disk/corona and single-phase ADAF, 
suddenly drops to the marginally stable orbit. In such 
case, when decreasing accretion rate attains the value 
of TOovap. branch, the discontinuous jump in the corona 
strength occurs from /cor < 1 to /cor = 1. This sharp 
transition is marked as a horizontal line called evapora- 
tion branch. In time-dependent solution the disc structure 
may actually follow this branch if m decreases suddenly. 

On the other hand, if we allow for existence of the 
reconstructed branch of the disk at the innermost part of 
the disk, whenever it is allowed, the solution enters the 
evaporation branch gradually, it is shallow and positioned 
at much lower accretion rate (long dashed line). However, 
such a solution has 'an ADAF hole' for a certain range of 
radii. In further considerations we adopt the first solution 
as the correct one since the inner secondary disk formation 
is not consistently described within the frame of our model 
and may not represent a physically acceptable solution. 

The mechanism of disk evaporation leads to departure 
from standard optically thick solution in the form of prac- 
tically horizontal branch which develops at intermediate 
accretion rates. For a — 0.1 the position of this branch is 
approximately at the same value of the accretion rate as 
the maximum of the accretion rate at the optically thin 
branch, as described by Narayan & Yi (1995b), Abramow- 
icz et al. (1995) and Chen et al. (1995), but this is not the 
case for other values of a since the optically thin branch 
was computed for a fixed assumed value of the Compton 
parameter characterizing the cooling of the optically thin 
flow. Correct solution should be obtained taking into ac- 
count the actual soft X-ray flux coming from outer parts 
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Fig. 16. The relation between the accretion rate and the 
surface density of the disk/corona system at 10i?schw for 
the viscosity parameter a = 0.1 in case of irreversible 
ADAF transition (continuous line) and with secondary 
disk rebuilding (long dashed line) . Short-dashed line shows 
the standard Shakura-Sunyaev model supplemented with 
advection essential at high accretion rates, calculated for 
M = IO^Mq - slight wiggles of the gas dominated lower 
branch are caused by bound-free opacities. Continuous 
line in the lower right part of the diagram shows solution 
for an optically thin flow (Narayan & Yi 1995b, Abramow- 
icz et al. 1995 ), computed using the code of Zdziarski 
(1998). 



of the disk and should be a natural extension of the two- 
phase solution in the case of cold disk absence. At present 
we cannot extend our horizontal branch to the low surface 
density range and to reproduce the solution for a single 
phase medium since the present model is based on as- 
sumption that the underlying disk is optically thick and 
absorbs half of the X-ray radiation emitted by the corona 
at a given radius and reemits this radiation in the form 
of soft flux. The description of the transition from a two- 
phase flow to a single-phase optically thin flow would re- 
quire relaxation of this assumption and a number of other 
modifications, like non-local description of the soft photon 
radiation field since the outermost part of the ADAF flow 
will be predominantly cooled by soft photons emitted by 
the outer parts (e.g. Esin et al. 1997) and/or soft photons 
originating from synchrotron radiation. We plan to study 
this problem in the future, in order to give a complete 
view of the transition between those two branches. 

The position of evaporation branch in our model de- 
pends on the viscosity - lower viscosity parameter a corre- 
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sponds to lower accretion rate at this branch. For a = 0.1 
and large black hole mass it joins the (almost) standard 
Shakura-Sunyaev branch just below the turning point 
marking the gas pressure dominated branch. Therefore, 
the transition practically happens from unstable radiation 
pressure dominated branch to evaporation branch with a 
very narrow range of accretion rates in between. Adopt- 
ing lower value of viscosity, however, broaden this range 
considerably. For low mass of the black hole the stable 
gas dominated branch exists even for a = 0.1 since the 
position of the transition from radiation pressure to gas 
pressure dominated solution in Shakura-Sunyaev disk de- 
pends on the mass of the c entral object. 

As seen from Fig. |l5|, the position of the evaporation 
branch does not depend on radius in the innermost part 
of the disk. At i? < i?ovac the position of the evaporation 
branch is given by Eq. p4| w hile for R > -Revac this posi- 
tion is determined by Eq. |2l|. This property has important 
consequences for the limit cycle evolution expected for a 
range of accretion rates corresponding to the unstable ra- 
diation pressure dominated branch of thick disk solution. 

3.3.2. Stationary solutions and outbursts 

This diagram (Fig. |l^, although not complete, shows the 
importance of the evaporation branch for stationary so- 
lutions and time-dependent evolution of accretion disk. 
If the accretion rate is larger than the Eddington limit 
the model is located on the stable optically thick ADAF 
branch. Stationary solution well represents an accretion 
process at this parameter range. When the accretion rate 
is very low, below the evaporation branch, the flow is de- 
scribed by the optically thin stable ADAF solution in the 
innermost part, and a disk/corona flow in the outer part. 
The solution should be also well represented by a station- 
ary model, unless the disk/corona interaction lead to some 
instabilities. 

However, for intermediate accretion rates, the 
disk/corona flow, as the classical radiation pressure dom- 
inated model, may display a complex limit cycle, as al- 
ready studied, without taking into account the evapora- 
tion branch, by Szuszkiewicz & Miller (1998), Honma, 
Matsumoto & Kato (1991) and Janiuk et al. (2000). The 
solution would oscillate between the upper branch during 
the eruption and and the evaporation branch between the 
outbursts. Actual position of the upper branch may be 
at lower accretion rates, if the strong outflow in addition 
to radial advection, is included (e.g. Nayakshin, Rappa- 
port & Meha 1999, Janiuk et al. 2000). The timescale 
for such evolution is shorter than the viscous timescale 
of the standard Shakura-Sunyaev stationary disk at this 
radius because of the reduced surface density. The exact 
prediction of the global outburst involves the entire unsta- 
ble part of the disk so, without actual computations, the 
exact timescale and the amplitude cannot be determined 



but order of magnitude timescales involved are seconds for 
galactic sources and months/years for AGN. 

The character of this evolution would depend signifi- 
cantly on the radial extension of the instability zone. If 
this extension is below i?evac the position of the evapora- 
tion branch is the same for all radii and the entire disk will 
oscillate between the upper (advection or outflow) domi- 
nated branch and the evaporation branch so the optically 
thick disk will be present during the entire cycle. However, 
if the radial extension is larger than i?ovac the evaporation 
branch at the outermost radius involved in the outburst is 
lower than evaporation branch inside. The accretion rate 
in low luminosity phase of the outburst is then too low for 
optically thick disk in the innermost part, the transition 
to ADAF takes place so the optically thick disk will be 
temporarily evacuated during such a cycle. It is interest- 
ing to note that probably both kinds of outbursts are seen 
in the case of microquasar GRS 1915-1-105, as two distinct 
kind of spectral behaviour and of the relation between the 
length of outburst and the frequency of QPO phenomenon 
are observed (Trudolyubov, Churazov & Gilfanov 1999). 
Also the observed X-ray delay in the outburst of soft X- 
ray transient GRO J1655-40 shows that most probably the 
system made the transition from an inner ADAF/outer 
disc/corona state to a pure disc/corona state Hameury et 
al. (1997). The same transition is seen in Cyg X-1 (Esin 
et al. 1997). 

The expected amplitude of the outburst is consider- 
ably reduced by the presence of the evaporation branch 
- expected luminosity changes are by an order of magni- 
tude while, if the gas pressure dominated branch is the 
lower branch of the limit cycle the outburst of up to four 
orders of magnitude would be expected, against any ob- 
servational evidence for the timescales involved. 

According to our model, there are strong systematic 
differences expected between the AGN and GBH since the 
range of the radiation pressure instability depends on the 
central mass while the position of the evaporation branch 
does not. If the evaporation branch is positioned signifi- 
cantly below the radiation pressure dominated branch we 
can also expect a stable accretion flow in the accretion 
rate range corresponding to the gas pressure dominated 
solution. It means that, independent on the viscosity, the 
range of accretion rates corresponding to stable, gas dom- 
inated solution is always broader for galactic black holes 
than for AGN. In particular, no such range at all is pre- 
dicted in AGN for a = 0.1 since this solution is entirely 
replaced by evaporation branch and ADAF. 

4. Discussion 

4-.1. Comparison to other models 
4.1.1. Underlying assumptions 

Our approach to modeling the disk/corona structure fol- 
lows the MMH94 attitude. We also model the energy dis- 
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sipation within a corona as due to the accretion flow, we 
take the transition between the disk and the corona which 
allows to determine the fraction of energy generated in the 
corona as a function of radius as well as the transition ra- 
dius to the optically thin flow as a function of the mass of 
the black hole, accretion rate and the viscosity parameter 
in the corona. 

There are, however, major important differences be- 
tween our model and MMH94: (i) our corona is basically 
a two-temperature medium with energy dissipation due to 
ion collisions and the Coulomb transfer of energy from ions 
to electrons, as in SLE and ADAF models (ii) apart from 
bremsstrahlung we allow also for Compton cooling (iii) we 
determine the conductive flux at the basis of the corona for 
the adopted heating/cooling mechanism instead of using 
the scaled prediction from solar corona. This last modifi- 
cation is responsible for the dependence of the accretion 
rate in the corona on the viscosity parameter in the case of 
our model and the lack of such dependence in MMH94 and 
subsequent papers (Liu, Meyer & Meyer-Hofmeister 1995, 
Liu et al. 1999, Meyer, Liu & Meyer-Hofmeister 2000). 
The predicted dcpcnd(nic;e of the coronal accretion rate 
on the radius is qualitatively similar in both papers in the 
case of low accretion rate but not in the case of large ac- 
cretion rate when the condensation, neglected by MMH94, 
is important. Our approach to conduction is therefore the 
same as of Dullemond (1999) but he did not consider a 
two-temperature plasma and Compton cooling. 

Two-temperature accreting corona was discussed by 
WCZ97 but their condition of disk/corona transition 
based on analogy with Compton heating medium was not 
accurate numerically and missed an important insight on 
corona formation in the innermost and outermost parts of 
the flow. 

4.1.2. Transition radius to ADAF-type flow 

The transition from an outer optically thick disk flow to an 
inner optically thin ADAF was predicted by a number of 
models. In some models the transition was based on ADAF 
principle (i.e. the transition happens whenever an ADAF 
solution is possible), but in other cases the necessity of 
the transition resulted from the model itself (MMH94 and 
applications, our model). It is therefore interesting to com- 
pare quantitatively those estimates as they may in princi- 
ple be estimated observationally. 

Condition based on the strong ADAF principle, i.e. on 
the absence of ADAF solutions: 

nr = 1.9 X 10%-2a^ liJschw (26) 

(Abramowicz et al. 1995, Honma 1996, Kato & Nakamura 
1998). 



Approximation to the results of Esin et al. (1997) for 
a = 0.25 and (3 = 0.5: 

= xl0*i?schw m < 0.084 

rtr = 30i?schw m « 0.084 ^^^^ 

= 3 - 30i?schw 0.084 < m < 0.092 
r-tr = 3iischw m > 0.092 

Liu et al. (1999) gives: 

rtr = 18.3m-i/^"i?schw (28) 
Our results: 

= 4.51m-i-33a5^i?schw m < 6.92 x IQ-^aH ^^^^ 

r-tr = 3iischw rh > 6.92 x lO-^a^ f 

All those conditions are roughly in agreement with ob- 
servational requirement that transition between the high 
and low state happens at luminosity 10'^^ erg s~^ (Tanaka 
1999). The accurate dependence of the transition radius on 
the accretion rate can in principle be determined from the 
shape of the reflected component (Done & Zycki 1999) in 
X-ray novae but the errors are considerable and the pres- 
ence of the hot Compton heated skin on the disk surface 
further complicates the interpretation of the data. How- 
ever, high precision spectroscopy from new X-ray satellites 
should give stronger constraints on the available models. 

4-2. Observational consequences 
4.2.1. Luminosity states 

The transition from optically thick flow to ADAF, i.e. the 
accretion rate at the evaporation branch marks the transi- 
tion between the hard state and soft state in accreting bi- 
naries, according to a number of suggestions (e.g. Narayan 
1996, Esin et al. 1997, Zdziarski, Lubihski & Smith 1999). 
When the accretion rate is high enough and two-phase 
flow exists up to the marginally stable orbit then the soft 
flux is always dominant, cooling of corona is efflcient and 
produces soft X-ray spectrum. Such a solution is stable 
in GBH for moderate accretion rates and may correspond 
to the soft (high) state. If the accretion rate is very high 
the disk is unstable and evolving between the evaporation 
branch and advection-dominated optically thick branch 
- such a parameter range may be appropriate for a Very 
High State, with its complex spectral behaviour (e.g. Tak- 
izawa et al. 1997, Revnivtsev et al. 2000). On the other 
hand, when the accretion rate is very low, the innermost 
part of the disk is an optically thin ADAF, dominating 
the spectrum, so such a case corresponds to the Low (or 
Hard) State. 

In our model, the position of horizontal evaporation 
branch does not depends on the mass of the black hole 
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revealed that in the intermediate state the reflection com- 
ponent is seen with an amphtude of order of 1 while the 
kinematical determination of the inner disk radius from 
Doppler shifts of the reflected component show disk dis- 
ruption (Zycki, Done & Smith 1998). We can interpret 
such a state within the frame of our model as correspond- 
ing to an accretion rate for which a strong corona forms 
at intermediate radii. In this case the reflection is strongly 
biased towards larger radii instead of towards innermost 
disk radius. Also very high values of the ionization pa- 
rameter may be consistent with the irradiating source be- 
ing directly above an accretion disk. Our interpretation 
is therefore different from that of Esin et al. (1997) since 
in their approach the intermediate state corresponds to a 
disrupted disk with an innermost ADAF flow instead of 
the two-phase disk/corona flow. 

4.2.3. Variability 



Fig. 17. The relation between the accretion rate and the 
surface density of the disk/corona system at 10i?schw for 
the viscosity parameter a = 0.1 in case of irreversible 
ADAF transition for an AGN {M — IOq, continuous line) 
and a GBH (M = 10©, dashed line). Stable solutions (pos- 
itive slope of the curve) are expected for intermediate ac- 
cretion rates in the case of a small mass of the black hole 
but this region is very narrow for massive black holes. 

(Fig. |l^). Nevertheless the turning point between the un- 
stable radiation pressure branch and the stable evapora- 
tion branch is different for AGNs and GBHs and the pos- 
sibility of formation of the High State in AGN is strongly 
suppressed. 

Since most of the energy generates at about 10 i?schw 
we can roughly identify the accretion rate ranges for var- 
ious luminosity states on the basis of Fig |l^. Assuming 
a = 0.1 wc obtain: 

State GBH AGN 

Very Hinh State to > 0.1 m > 0.07 

^ ^ (30) 

High State 0.07 < to < 0.1 

Low State to < 0.07 to < 0.07 

However, we stress that those limits depend signifi- 
cantly on the adopted value of the viscosity parameter a. 

4.2.2. The nature of the Intermediate State in GBH 

The transition of X-ray novae from soft state with the 
dominant disk emission to a hard state with the dominant 
thin plasma (e.g. ADAF) contribution proceeds through 
the Intermediate State. The analysis of the Nova Musca 



The lack of High State in AGN, as predicted by our model 
for the viscosity parameter a = 0.1 may offer an explana- 
tion why the observed variability properties of AGN and 
GBH are not the same. GBH, when in their High State, 
dominated by disk emission, are only weakly variable in 
X-rays in the timescales of about a second, with rms am- 
plitude of order of a percent (e.g. Takizawa et al. 1997). 
Spectrally, such states are basically similar to quasar spec- 
tra. However, quasars are strongly variable, with rms am- 
plitude in the optical band dominated by the disk emis- 
sion is of order of 15 per cent (e.g. Giveon et al. 1999) 
in the timescales of years, so the simple scaling with the 
mass of the black hole does not seem to hold. It may mean 
that quasars are actually counterparts of the GBH in their 
Very High State. Since the spectra of GBH in Very High 
State are quite complex and not always dominated by the 
disk emission we can also expect such weak disk spectral 
states among the quasar/Seyfert class. Well studies source 
MCG-6-15-30 (Tanaka et al. 1995 and subsequent papers) 
may be an example, due to its exceptionally low mass and 
consequently high luminosity to the Eddington luminosity 
ratio (Nowak & Chiang 1999). 

4.2.4. Broad Line Region and LINERs 

Nicastro (2000) proposed a very interesting model of the 
Broad Line Region based on the disk/corona model of 
WCZ97. The basic idea is that the strong outflow develops 
at radii where the corona strength reaches its maximum 
and cool clouds forming in this disk wind are responsible 
for the emission of broad lines. Predictions based on our 
new model does not alter presented estimates significantly. 
The radius of the maximum corona strength in the new 
model can be determined from the condition to^ = (see 
Eq. |2^), using the result with negligible evaporation to 
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determine the local coronal accretion rate (Eq. which 
gives the dependence 



1.02 X lO'^m 



3^16/23 37/115 



0.1 



(31) 



smaller by a factor of few than the WCZ97 formula (see 
Eq. 3 of Nicastro 2000). On the other hand the new corona 
efficiency is now roughly symmetric around this point so 
the weighted formula for the outflow used by Nicastro 
would not decrease this radius, as it did in the case of 
sharp edge corona of WCZ97 so the net effect will be 
roughly the same as before. 

It was suggested by Lasota et al. (1996), that LINERs 
may represent the systems with ADAFs. The transition 
from objects with broad emission lines to LINERs hap- 
pens, according to Nicastro (2000), due to the switch off 
the outflow if the accretion rate is too low for the exis- 
tence of the radiation pressure dominated zone. The most 
efficient part of the corona, however, is located in the gas 
pressure dominated region so this connection is not clear. 
On the other hand, for very small accretion rate, within 
the frame of our new model, there is a transition from 
secondary disk rebuilding phase to a definite ADAF flow 
in the innermost part (see Eq. p3|) . Such a change may 
be important from the point of view of the radiative cool- 
ing, the efficiency of advection and, consequently, the mass 
outflow. The mass loss through the wind as a function of 
radius is not calculated yet in the present model but this 
can be done in the future, following the method of e.g. 
Meyer & Hoffmeistcr (1994) and WCZ97. 
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Appendix A: Evaporation/condensation rate at 
the basis of the corona 

The electron temperature profile within the transition 
zone is determined by the effect of heating, radiative cool- 
ing, conduction and the vertical motion of the material 
(e.g. Begelman & McKee 1990, McKee & Begelman 1990, 
MMH94). The general form of this equation in a plane 
parallel geometry is the following 



d_ 

dz 



5 fcTi 
2 mn 



dF, 



cond 



dz 



Qt 



A, 



(A.l) 



where the conduction flux due to electrons is given by 

(A.2) 



^ cond — f^Q-^c T 

dz 



Here we included the effect of radial advection into the 
cooling term and we neglected the conduction due to ions. 



Following McKee & Begelman (1990) we multiply the 
Equation A.l by fcond and integrate in the vertical direc- 



tion across the transition zone: 
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5 kTi ' 


Jzi dz 


2 mn . 



2 [-P'cond(^2) 



Te(P) 



[Qt - K]n,T^J^dT, 



cond(^l)] — 

(A.3) 



The conduction flux at the upper and lower point of 
the transition zone vanishes which is used as a bound- 
ary condition in full numerical computations (e.g. DuUc- 
mond 1999) so the second term on the left hand side of 
the equation vanishes thus leaving a relation between the 
mass transfer (flrst term) and the integral on the right 
hand side. 

We assume that the transition from hot coronal to cold 
disk material in the vertical direction takes place in a very 
narrow geometrical zone so any change of the pressure 
across this zone can be neglected. This allows to calculate 
the right hand side integral conveniently, using algebraic 
expressions for the density, ion temperature and the heat- 
ing and cooling as functions of the electron temperature. 

The left hand side can be approximately expressed by 
estimating first the Field length (Field 1965) measuring 
the thickness of the transition zone 



Af = 



<\Qt-A>J 



|/^"^'''[Qe+-A]KoTe'/'dre|l/2' 



(A.4) 



We now estimate the typical value of the conductive flux 
within the transition zone as 

< Fcond >- KoT,iPy/yXF (A.5) 

which allows us to express the left hand side term as 

5 kTdP) 

^nddz Fcoiid > m^-z — -(A. 6) 

2 mn 

Therefore the approximate formula for the condensa- 
tion/evaporation rate at the disk/corona boundary reads 
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5kTi' 
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